Zooming into the coexisting regime of ferromagnetism and superconductivity in 

ErRh4B4 single crystals 
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High resolution measurements of the dynamic magnetic susceptibility are reported for ferro- 
magnetic re-entrant superconductor, ErRh4B4. Detailed investigation of the coexisting regime 
reveals unusual temperature-asymmetric and magnetically anisotropic behavior. The supercon- 
ducting phase appears via a series of discontinuous steps upon warming from the ferromagnetic 
normal phase, whereas the ferromagnetic phase develops via a gradual transition. A model based 
on local field inhomogeneity is proposed to explain the observations. 

PACS numbers: 74.25.Dw; 75.50.Cc; 74.25. Ha; 74.25.-q; 74.90.-|-n 



The coexistence of the long-range magnetic order and 
superconductivity was first discussed even before the ap- 

Sarance of the microscopic theory of superconductivity 
. Since then this topic remains one of the most inter- 
esting and controversial in the physics of superconduc- 
tors with many reviews and books devoted to the subject 
0, [3, 0, S 0| • Despite significant effort in new materi- 
als design and discovery, there are only few, confirmed, 
ferromagnetic superconductors. Local, full-moment fer- 
romagnetic superconductors: ErRh4B4 (Tfm ~ 0.9 
K, « 8.7K), Ho.MoeSg H {TpM ~ 0.7 K, « 1.8 
K); the weakly ferromagnetic ErNisBaC 0] {TpM ~ 2.3 
K, Tapm « 6 K, Tc ~ 11 K); and more recent itin- 
erant superconducting ferromagnets, UGe2 {TpM ~ 33 
K, Tc « 0.95 K at P w 1.3 GPa) and URhGe 
{Tfm ~ 9.5 K, Tc « 0.27 K). Whereas all these ma- 
terials are very interesting on their own, the coexistence 
of growing, large, local moment ferromagnetism and su- 
perconductivity is most clearly presented in ErRh4B4, 
which is the subject of the present Letter. In particular, 
we are interested in the details of the narrow tempera- 
ture interval (~ 0.3 K) where the two phases coexist and 
infiuence each other. 

ErRh4B4 was extensively studied over past 30 years 
0) H, 0, S 0- The ferromagnetic phase is primitive 
tetragonal with the c— axis being the hard and a— axis 
being the easy magnetic axes. Detailed measurements of 
anisotropic magnetization and upper critical filed, 11^2 
were done by Crabtree et al. [ll|, [l3| , who found that 
(along a— axis) peaks at 5.5 K due to large param- 
agnetic spin susceptibility in that direction [1, 11 1 . In the 
contrast, collapses near the onset of the long-range 
ferromagnetic order. 

Neutron diffraction studies have established the ex- 
istence of a modulated ferromagnetic structure at the 
lengthscale of ~ 10 nm [l3, 14|. In single crystals, results 



suggested that coexisting phases consists of a mosaic of 
normal FM domains and SC regions larger than ^ 200 
nm in size. The SC regions contain modulated FM mo- 
ment with a period of ~ 10 nm. These regions could be 
regular domains, spontaneous vortex lattices or laminar 



structures with > 200 nm periodicity and modulated SC 
domains in between [3]. Thermal hysteresis is observed 
both in the normal Bragg peak intensity and the small- 
angle peaks. For the small-angle peaks, the intensity is 
higher on cooling than on warming. This is opposite to 
the behavior of the regular Bragg peaks from the FM 
regions [l^. Furthermore, the first-order transition, ob- 
served in satellite peaks temperature dependence |l4l|. i s 
consistent with the spiral state of Blount and Varma [l5| . 
However, a continuos transition was reported in other 
neutron diffraction l^, 1^ and specific heat experiments 
[13] ■ Such a transition can be realized in a modulated 
structure or via spontaneous vortex phase. 

Theoretically, some striking features of the coexist- 
ing phase include an inhomogeneous, spiral, FM struc- 
ture jlSl. Il8| or a fine domain, " cryptoferromagnetic" 
phase [2, ll9|, a vortex - lattice modulated spin struc- 
ture [2o| . type-I superconductivity 0, [IS HH , a gapless 
regime and possibly, an inhomogeneous Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) state 0]- Another interest- 
ing possibility is the development of superconductivity at 
the ferromagnetic domain walls [l^, [31 ■ 

In this Letter we report precision measurements of the 
dynamic magnetic susceptibility of ErRh4B4 with an em- 
phasis on the narrow temperature region where ferromag- 
netism and superconductivity coexist. We find that the 
transition is highly asymmetric when FM — + SC (heat- 
ing) and SC FM (cooling) data are compared. The 
FM^SC transition proceeds via a series of discrete steps 
from FM to SC phase upon warming and proceeds via a 
smooth crossover from the SC to FM state upon cooling. 
With this new information we analyze relevance of some 
predictions made over years for the coexisting phase. 

Single crystals of ErRh4B4 were grown at high temper- 
atures from molten copper fiux as described in [3, [11] ■ 
Resulting samples were needle shaped with crystallo- 
graphic c-axis along the needle. Transport measurements 
gave residual resistivity ratio of about 8, consistent with 
previous reports. The anisotropic Hc2 (T) curves (see in- 
set to Figj3] below) are consistent with earlier reports as 
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FIG. 1: Dynamic magnetic susceptibility in ErRh4B4 single 
crystal measured along the magnetic easy axis (perpendicu- 
lar to the needle-shaped sample). Each curve corresponds to 
a fixed value of the applied dc magnetic field in the range 
indicated in the figure, (color online) 
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FIG. 2: Magnetic susceptibility in a ferromagnet- 
superconductor transition region measured at different ap- 
plied fields. Note the temperature scale and highly assymet- 
ric character of the FM-SC and SC-FM transitions, (color 
online). Blue - warming, red - cooling. 



The AC magnetic susceptibility, %, was measured with 
a tunnel-diode resonator (TDR) which is sensitive to 
changes in susceptibility Ax ^ 10~^. Details of the mea- 
surement technique are described elsewhere 
In brief, properly biased tunnel diode compensates for 
losses in the tank circuit, so it is self-resonating on its 
resonant frequency, uj = l/VXC ^ 10 MHz. A sample 
is inserted into the coil on a sapphire rod. The effective 
inductance changes and this causes a change in the res- 
onant frequency. This frequency shift is the measured 
quantity and it is proportional to the sample dynamic 
magnetic susceptibility, x 21 > 11] ■ Knowing geomet- 
rical calibration factors of our circuit, we obtain x (T, H). 
Advantages of this technique are: very small AC excita- 
tion field amplitude (~ 20 mOe), which means that it 
only probes, but does not disturb the superconducting 
state; high stability and excellent temperature resolution 
(~ 1 mK), which allowed detailed study of the coexist- 
ing region, which is only ~ 500 mK wide. Normal-state 
skin depth is larger than the sample size, so we probe the 
entire bulk in the coexisting region, but when supercon- 
ducting phase becomes dominant, there is a possibility 
that some FM patches still exist, but are screened. 

Figure[T]shows the full temperature scale magnetic sus- 
ceptibility, X, in single crystal ErRh4B4 for an applied 
field oriented along the easy axis and perpendicular to 
the needle-shaped crystal. The peak in x {T) at the fer- 
romagnetic to superconducting boundary below 1 K is 



a common feature observed in local moment ferromag- 
nets (29| . Clearly, superconductivity is fully suppressed 
in the ferromagnetic phase. Note that at elevated fields, 
the response is nonmonotonic on the SC side close, to 
the FM boundary, indicative of enhanced diamagnetism 
(larger, negative x), which may be due to suppressed 
magnetic pairbreaking or entering into another phase, 
such as FFLO 3- 

Figure [2] zooms into the SC^FM transition region. 
Measurements were taken after zero-field cooling, apply- 
ing external field and warming up (ZFC-W) above Tc 
and then cooling back to the lowest temperature (FC-C). 
There is striking asymmetry of the transitions - when the 
superconducting phase develops out of the FM state, the 
response proceeds with jumps in the susceptibility, which 
are clearly associated with the appearance of supercon- 
ducting regions of finite size. The steps are present up 
to the largest field at which superconductivity survives. 
Decreasing temperature shows a completely different re- 
sult: the transition is smooth and gradual and proceeds 
to lower temperatures. 

To better understand the dynamics of the transition, 
the top frame of Fig.[3]shows measurements at H — for 
different temperature ramp rates. Temperature variation 
is shown in the inset. These data clearly demonstrate 
that this hysteresis is insensitive to heating/cooling rates. 
It should be noted that all the other data presented in 
this Letter were taken with the slowest cooling rate of 60 
/iK/s. 

Similar hysteresis and steps are also present in the an- 
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FIG. 3: Hysteresis of the transition at zero applied field mea- 
sured at diflterent ramp rates. Temperature sweep profiles are 
shown in the inset. Clearly, this hysteresis and jumps are not 
dynamic effects. 
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FIG. 5: Details of the hysteresis with partial scans as ex- 
plained in the text. Inset shows anothe partial loop on a 
cooling part of the curve. 
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FIG. 4: Transition region measured for magnetic field ori- 
ented along the needle and c-axis. Inset: summary phase di- 
agram for two orientations measured by transport (symbols) 
and TDR. 



other orientation, when external magnetic field is applied 
along the c— axis. This is shown in Fig. 21 Note that peak 
in X (T) at the FM boundary is not present, which is 
consistent with the behavior of anisotropic local-moment 
ferromagnet [2^ . The inset to Fig. [4] shows the phase di- 
agram obtained from resistivity and TDR measurements 
for both orientations. There is excellent agreement be- 
tween the two techniques and, as noted earlier, this dia- 



gram is consistent with previous reports [ll|, ■ 

Finally, Fig. [5] shows so called minor hysteresis loops 
(not as function of field, but temperature). The labels 
show the evolution of the susceptibility. It starts from 
low temperature at (1) when sample was warmed up to 
first signs of superconductivity that appeared as small 
jump at (2), then warmed further and reaching almost 
full superconductivity at (3), but then cooled back down 
to (4) as indicated by arrow and warmed back to (5). 
Note that along (3) (4) x is significantly different 
from (4) (5). Another similar minor cooling- warming 
loop follows (6) (7) (8) after which the sample was 
cooled down to return to (10) = (1) via (9). Interestingly, 
there are no steps or jumps observed on the minor loops 
even on warming. Also, the slope dx/dt is similar on cool- 
ing and warming and is very different from the original 
steep slope (2) (3). This is consistent with the present 
of vortices, probably pinned by the modulated FM/SC 
structure. The inset in Fig. [5] shows a small minor loop 
on a cooling part. This loop has small slope comparable 
to the larger loops described in the main frame. 

Let us now turn to the interpretation of these results. 
Clearly, the FM^SC transition proceeds via a series of 
jumps in diamagnetic screening due to formation of su- 
perconducting regions of macroscopic volume, roughly 
5% — 20% of the sample volume depending on the ap- 
plied field and temperature. Indeed, each observed step 
may be a result of simultaneous formation of many in- 
dividual superconducting domains of similar size. These 
steps in x are present both for H\\c and H\\a axes, al- 
though in the latter case the steps are smaller and are 
more pronounced, possibly due to magnetic and shape 
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anisotropies. The number of steps increases with the in- 
creasing field and the first step (the first sign of supercon- 
ductivity) occurs at a higher temperature for larger ap- 
plied field. Overall, FM^SC transition is apparently of 
the first order and exhibits behavior consistent with t ype - 
I superconductivity as predicted theoretically 0, 2^ 21 1 . 

From our point of view, the first jump occurs at a tem- 
perature where internal field is equal to a supercooling 
field of a type-I superconductor. (Note that apparent su- 
perheating of the FM^SC transition corresponds to 
supercooling of a regular type-1 superconductor 



17, 3 



because we enter normal phase on cooling). When first 
superconducting domains appear, effective magnetic field 
around them increases due to the flux expulsion and in- 
ternal field becomes more inhomogeneous. This net in- 
crease in the internal field in the remaining FM regions 
stabilizes them to higher temperatures. The system now 
needs to get farther away from the initial FM boundary, 
deeper into the SC state to produce more superconduct- 
ing patches. In this scenario, the observed jumps in x cor- 
respond to a cascade of supercooling transitions. If the 
temperature is lowered before the transition is complete 
domains stay stable to lower temperatures due to physics 
similar to superheating of a type-I superconductor. It is 
also quite possible that superconducting domains have 
the modulated spin structure seen in neutron scattering 
Q. Finally, it seems that ferromagnetic domains are not 
directly related to the observed steps, because at higher 
fields, the number of these domains decrease and domi- 
nant domains (along the applied field) grow in size. 

In a striking contrast with FM^SC transition, the 
SC— >FM transition is smooth and proceeds to much lower 
temperatures. Yet the transition is hysteretic as evident 
from the minor loops shown in Fig. [5l It is possible 
that Abrikosov vortices are being spontaneously created 
as the temperature is lowered and the systems crosses 
over into the normal state when vortex cores overlap. 
The vortex state is also compatible with long-range co- 
herence observed in neutron scattering experiments [l^ . 
At the same time ferromagnetic modulation with a pe- 
riod of ~ 10 nm may also develop between the vortices 
0, [H, [H, [T^. This would be also be similar to FFLO 
state in the presence of vortices [sij. Moreover, this 
would explain different intensities of small-angle satel- 
lite peaks, because coherence volume in the vortex state 
must be much larger compared to domain-like state on 
warming. We also note that unusual enhancement of dia- 
magnetism in the vicinity of the FM boundary from the 
SC side, could be due to an FFLO pocket as predicted by 
Bulaevskii for ErRh4B4 ■ If we plot temperature of a 
minimum in x (T) as function applied field, and also i?c2, 
we obtain phase diagram remarkably similar to Fig. 7 of 
Ref. 2]. 

Overall, we conclude that we are witnessing an unusual 
transition. It is definitely first order on warming with 
signs of type-I superconductor "supercooling". However, 



it is smooth upon cooling and exhibits smooth minor 
loops similar to a type-II superconductor. A second or- 
der transition occurs between normal and FFLO phases 
as well as between normal and SC for type-II supercon- 
ductor. However, transition from SC to FFLO state is 
first order as well as from SC to spiral state. It is possi- 
ble that size of the new phase nuclei is so small that we 
cannot resolve it upon cooling. 

Discussions with Lev Bulaevskii, Alexander Buzdin, 
Vladimir Kogan and Roman Mints are appreciated. 
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